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Abstract 

We propose a novel mechanism for reproducing the reahstic hierarchical structure of the 
observed CKM mixing matrix and quark masses by means of introducing a warped metric. We 
illustrate the method on the basis of a specific Type IIA intersecting brane model proposed 
by Ibafiez et al. The model is compactified on a direct product of three tori, on which four 
D-brane stacks wrap, and the intersections between two stacks collectively provide three gen- 
erations of quarks and leptons. The exponential of the area formed by the intersection points 
corresponding to three matter fields is interpreted as the Yukawa coupling among them. It is 
known, however, that the Yukawa matrices for quarks are generically singular, and that it does 
not give definite mixing angles nor hierarchically suppressed quark masses. We show that the 
newly introduced warp effect on the internal manifold modifies the Yukawa matrix elements and 
generate hierarchical quark masses and CKM mixing angles. 
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1 Introduction 



The string theories have been regarded as one of the most promising models for describing both 
gravity and particle physics. Many researchers have tried to derive the Standard Model within the 
framework. Since the theories are originally constructed in higher dimensional spacetime, the extra 
dimensions must be compactified on some manifold. Through these procedure, in most cases, there 
appear many scalar fields at low energy scale, and their vacuum expectation values remain unfixed 
because we cannot calculate the potentials of these fields within the perturbative framework. 

In the mid-1990s, the author in Ref. [1] formulated D-branes, spatially extended objects, on 
the worldsheet of string theories by imposing Dirichlet boundary conditions on the end points of 
strings. It is found that D-branes carry a certain amount of Ramond-Ramond charge, and that 
they could be regarded as solitons in supergravity theories. The emergence of solitons, as is known 
from quantum field theories, suggests a duality transformation from one theory described in terms 
of a certain set of fundamental fields and couplings to another in different terms. This string duality 
shed light on the relation among string theories, and invoked a considerable number of studies on 
the construction of phenomenological models with D-branes in Type IIA, Type IIB, and Type I 
string theories although they are once regarded inadequate for explaining the real world. 

Recently, intersecting D-brane models have been intensively studied. The scenario was first 
proposed by Blumenhagen et al. in the context of Type I string theory [2, 3, 4]. On a stack of 
N D-branes, the charges assigned to the end points of open strings, called Chan-Paton factors, 
are known to form collectively the adjoint representation of U{N) gauge group. Furthermore, at 
the intersection between a stack of N D-branes and that of M D-branes there appears a massless 
fermionic field belonging to a bifundamental representation {N, M) of a gauge group U (N) x U (M) 
[5, 6] . In the framework of intersecting brane models, these massless open string modes are identified 
with matter fields such as quarks and Icptons. It is worth noting that a Yukawa coupling is described 
as e"*^"^, where /c is a constant of order 1, and A is the area formed by three intersecting points 
corresponding to the relevant matter fields. Since the superstring theories are constructed in 10 
dimensional spacetime, we have to compactify them on 6 dimensional manifold when it comes to 
discussing phenomenology. In that case, D-branes wrap on some cycles of the manifold and they 
generically intersect with each other several times, and the intersection numbers determine the 
number of matter fields. Many researchers have investigated intersecting brane models in various 
aspects: the search for the standard model [7, 8, 9, 10, 11] , the construction of grand unified 
theories[12, 13, 14, 15], supersymmetric models [16, 17, 18, 19] and composite models[20, 21], the 
CFT calculation of couplings between fields locating at intersections [23, 24, 25], quark masses and 
a CKM mixing matrix [11, 21, 22, 26], and phenomenological viability [27, 28, 29]. 

In this paper, we focus on a specific Type IIA intersecting brane model compactified on a direct 
product of three tori proposed by Ibahez et al.[9] In the model, the family of right-handed quarks 
is distinguished by the difference of intersecting points in position on the first torus, and that of 
left-handed ones on the third torus. The Yukawa matrices are described by the area formed by 
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three intersecting points corresponding to quarks and a Higgs boson. Owing to the direct product 
structure of the internal manifold, the Yukawa matrices take the form of CL-ii^^^ and have only 
one non-zero eigenvalue, which is identified to be the third generation quark mass. In addition, 
we cannot obtain a definite CKM mixing matrix in the model. In this paper, we investigate 
quark masses and a CKM mixing matrix assuming a warp effect on the internal manifold, i.e., the 
dependence of the first torus volume on a coordinate of the third one. We reproduce the realistic 
hierarchical structure of the observed CKM mixing matrix and quark masses under some condition 
on small parameters originated from the warp effect and the size of tori. 

The outline of our paper is as follows. In section 2, we briefly review the construction of 
intersecting brane models. In section 3, assuming a warp effect on tori we calculate the Yukawa 
matrices and obtain the mass spectra for quarks and the CKM matrix. We show that our model 
reproduces the realistic hierarchical structure with several small parameters determined properly. 
The section 4 is devoted to the discussion. 

2 Intersecting Brane Models 

We briefly review an intersecting brane model presented in Ref . [9] . We consider Type IIA string 
theory compactified on orientifold six tori, with several stacks of D-branes intersecting each other 
at angles. Let us assume that orientifolds are located at = 0, tt, = 0, tt and xs = 0, tt, and that 
D6-branes wrapped on a direct product of 1-cycles of each torus with a set of winding numbers 
(n^,m^). The intersection is four dimensional. We know that at the intersection between Na D6- 
branes and Nj, D6-branes, there is an {Na,Nb) representation of the gauge group U{Na) x U{Nh). 
In our setup, there are also mirror branes owing to the orientifolds. The winding number of mirror 
branes are described as (n^,— m^). Later on, two sectors play an important role in generating 
quarks and leptons. The first sector Da-Dj,, realized on the intersection between two D-brane 
stacks, D6a-branes and DGft-branes. The intersection number turns out to be 

= {nlml-nlml){nlrnl-nlml){nlml-nlml), (2.1) 

and the sector gives a bifundamcntal representation (Na, iV^). The second sector Da-^D^ is realized 
on the intersection between D6a-branes and a mirror stack corresponding to another D-brane stack, 
Z)6fe-branes. The intersection number is 

lab = -inlml+nlml){nlml + nlml){nlml + nlml), (2.2) 

and it gives a representation {Na,Nj,). A negative sign in (2.1) and (2.2) represents the opposite 
chirality. The other sectors are not relevant to our discussion. 

So far, we have not imposed any restriction on winding numbers. We have to take account of 
Ramond-Ramond tadpole cancellation condition, which means the cancellation of D-brane charge 
on the internal manifold. If there is no orientifold, the condition is described as 

K 

llNana = 0, (2.3) 

a 
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i.e. the sum of the cycles on which D-branes wrap is equal to zero. If there are several orientifolds, 
they play a role of negatively charged objects. In our setup, the orientifold planes wrapping on a 
cycle (1,0), (1,0), (1,0), the condition reads 

^Nanlnlnl = 16, 

a 

Y^Namlmlnl = 0, 

a 

J2Namlnlml = 0, 

a 

Y^Nanlmlml = 0. (2.4) 

a 

We should consider the possibility of NS B-flux 6* on a torus T* [4] . Starting with Type I string 
theory compactified on a direct product of three tori, the introduction of NS B-flux background 
6 = 1/2 on a torus modifies its complex structure. In the intersecting brane picture, the effect 
is equivalent to changing a winding number on an original torus (n,m') into {n,m' + n/2). We 
will describe the configuration of D-branes in terms of this newly defined effective winding number 
m = m' -|- n/2 on a torus with non-zero S-field background. * 

We are now ready to look for winding numbers generating standard model matter fields: 

gi = (3,2,l/6); C/ij = (3, 1,-2/3); = (3, 1, 1/3); 

L^ = (1,2,-1/2); S|j = (1,1,1). 

Note that we have to put at least four stacks of D-branes. This is because the right-handed 
charged lepton Ej^ does not have the charge of SU (3) nor that of SU (2) , thus at least two more 
D-brane stacks are required. Prom now on, we discuss the models with four D-brane stacks: stack 
a composed of 3 D6-branes, stack 6 of 2 D6-branes, and the others composed of a D6-brane. 

First, the requirement of the above matter fields allows only the winding numbers denoted in 
Table 2, where /3* = 1 — 6*, e = ±1 and p = 1, 1/3. The variables n^, n^, and take integers. 
Note that on the third torus, the B field always takes a value = 1/2. This model contains four 
Z7(l) charges originated from D-brane stacks Qa, Qb, Qc, and Qd- Each of them can be understood 
in terms of the Standard Model framework: 

Qa = SB, Qc = 2Ir, Qd = -L, (2.5) 

where B is the baryon number, L is the lepton number, and Ir is the third component of right- 
handed weak isospins. Qi, plays a role of a generation dependent Peccei-Quinn symmetry. Owing 
to the couplings to the Ramond-Ramond fields, three of the four U{1) charges become massive and 
only the gauge boson corresponding to the hypercharge 

Qy = ^Qa - \Qc + \Qd (2.6) 

*The existence of orientifolds on tori requires the quantization of the B flux 6 = 0, 1/2. 
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remains massless if and only if the following relation holds: 



(2.7) 



Ni 


inj,mj) 


(nlmf) 




iVa = 3 


(l//3i,0) 




(l/p,l/2) 


N, = 2 


(ni,-e/3i) 


(l//32,0) 


(l,3p/2) 


N, = l 


(ni,3pe/3i) 


(l//3^0) 


(0,1) 


Nd = l 




{nl -(5^1 p) 


(l,3p/2) 



Table 1: The winding numbers yielding the Standard Model matter fields. 

Secondly, as for the tadpole cancellation, these winding numbers satisfy automatically all but 
the first condition of Eqs. (2.4). It reads 

16. (2.8) 



If there are more than four D-brane stacks, their contribution should be added to the left handed 
side of Eq. (2.8). It is possible to satisfy tadpole cancellation conditions by adding some D-brane 
stacks which do not intersect with the first four stacks, i.e. keeping the intersection numbers the 
same. 

In conclusion, the desired configuration is realized by the winding numbers defined in Table 2 
with the parameters satisfying Eqs. (2.7) and (2.8). The resultant matter contents and their U{1) 
charges are shown in Table 3. We are now ready to consider the Yukawa couplings. For simplicity, 
we take a specific D-brane configuration and calculate the Yukawa matrices 

= = 1, (2.9) 
e = -l, p=l, (2.10) 



2, 4 = 0, nl 



1, nl 



0. 



(2.11) 



The brane configuration of this setup is depicted in Fig. 1, where we take specific values for 
moduli parameters indicating positions of branes. In order to satisfy the first tadpole cancellation 
condition, we add 5 parallel D-branes with a winding number (1,0) (1,0) (1,0), which do not 
intersect with any D-brane stacks. 

In this context, the Higgs scalar fields come from the open string modes stretching from two 
parallel D-brane stacks not from an intersection. These modes are tachyonic, and are expected to 
undergo tachyon condensation, which could be liken to the Higgs mechanism, and to become finally 
ordinary Higgs fields with appropriate masses. Nevertheless, the mechanism of condensation has 
not been clarified. For the purpose of avoiding the complexity on the tachyon condensation, we 
assume that the size of the second torus is so small compared with the others that we neglect it 
when evaluating Yukawa couplings. 
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{n},mj) 




{nf, mf) 


Na = 3 


(1,0) 


(2,-1) 


(1,1/2) 


Nb = 2 


(0,1) 


(1,0) 


(1,3/2) 


Nc = l 


(1,-3) 


(1,0) 


(0,1) 


Nd = l 


(1,0) 


(0,1) 


(1,3/2) 



Table 2: A combination of winding numbers yielding the SM matter fields 



Intersection 


Matter fields 


repr. 


Qa 


Qb 


Qc 


Qd 


Y 


(a, 6) 


Ql 


(3,2) 


1 


-1 








1/6 


ia,b*) 


QL 


2(3,2) 


1 


1 








1/6 


(a,c) 


Ur 


3(3,1) 


-1 





1 





-2/3 


(a,c*) 


Dr 


3(3,1) 


-1 





-1 





1/3 


{b,d*) 


L 


3(1,2) 





-1 





-1 


-1/2 


icd) 


Er 


3(1,1) 








-1 


1 


1 


icd*) 


Nr 


3(1,1) 








1 


1 






Table 3: Matter contents and U{1) charges 

This model has the following Yukawa interactions: 

yfQLUih+yfQLDi^H2 
+ yr,qlUi,H,+yf^qiD^j,h2 

+ yijL'Ej^H2 + yfjVNjih,+h.c. (2.12) 

where Ql and q\^{i = 1, 2) are left-handed quarks, Uj^{j = 1, 2, 3) are right-handed ones, and hi, h2, 
Hi, and H2 are Higgs fields. The coefficients y^ , y^,y^, y'^, y^ , and y^ are Yukawa couplings. Each 
of them is described by the triangular area Aijj^ formed by three intersecting points corresponding 
to the fields included in the Yukawa interaction term 

yi3k = e-^^'i>', (2.13) 

where A; is a constant of order 1, which is expected to be fixed through stringy calculation [23, 24]. 

This expression can be interpreted intuitively as follows. All the matter fields such as quarks 
and leptons come from open string modes, with one end on a stack and the other end on the other 
stack. For example, Q\ has one end on D6a-branes and the other on DB^-branes. Likewise, U}^ is the 
open string mode stretching between DG^-branes and D6c-branes, and hi between D6c-branes and 
D6a-branes. The Yukawa coupling among these three fields evaluated by the scattering amplitude 
among these open strings, i.e., by the exponential of the area swept by them. When the compact 
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manifold is a direct product of tori, Yukawa couplings are expressed by the sum of the triangular 
areas projected on each torus j4j-"^[11]: 

A^,k = j:A'^l (2.14) 

n=l 

Before we come to the main subject, we must draw attention to the stability of intersecting 
brane models. It is known that, in general, parallel D-branes do not attract each other because 
of the cancellation between the repulsive force coming from RR sectors and the attractive one 
from NS sectors. However, as for D-branes intersecting at angles, this is not the case. Besides the 
massless fermion states, in general, there are some scalar states at each intersection, which could 
be regarded as supersymmetric partners of the fermion states. The mass square of these scalar 
states is parameterized by the angles between two D-brane stacks. Some choices of angles give a 
negative mass square, which suggest the instability of the configuration. However, with the winding 
numbers adopted here, we could choose the radii of tori with which tachyons do not appear at any 
intersection. 




Figure 1: The configuration of D-brane stacks on the first, second, and third torus (from left to 
right). The stack a is composed of 3 D6-branes, the stack 6 of 2 D6-branes, and the stack c is 
a single D6-brane. The broken lines and the dotted ones represent mirror branes and orientifold 
planes, respectively. Since DGc^-branes are not used for our analysis on quark sectors, they are not 
depicted. On the second torus, the open strings stretching between b-c* and b-c correspond to the 
i/j's, and the /i^'s. On the third torus, the unshaded portion is a fundamental region. 



3 The intersecting brane model with a warp effect 

We now show that a warp effect plays an important role in generating hierarchically suppressed 

(n) 

fermion masses and CKM mixings. First, let us consider the up-quark sector. We define Ajj -- as 
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the area of a triangle formed by three intersection points Qu, Unj and Xu^^ on the n-th torus, 
where Qu = Qi (i = 1, 2), Qls = Ql, Xjj,, = hi {i = 1, 2) and Xu^ = Hi. 

Then the Yukawa matrix for the up sector is described as = {Xij^i)e^^^^^-^^ , where Aij^ij = 
■^u\j + ^u\j- there is no warping, owing to the degeneracy of the locations of quark doublets 
on the first torus, depends only on the index j, where let us call it . On the other hand, 

on the third torus the right handed quarks stay at the same point, thus we define 
Therefore, the Yukawa matrix for up-sector quarks is described as 

Yu,ij = {Xu,i)aibJ (3.15) 

where Oj = e~^^? and b'j = e~^^^ . For the down-quark sector, we obtain 

YD,^j = {XD,i)aibf (3.16) 

where bj' = e~^^3 , and A^ is the area formed by Qli-, Drj and XD,i on the first torus, where 
XD,i = h2 {i = 1, 2) and Xd,3 = H2. For simplicity, we assume all the vacuum expectation values 
of Higgs fields to be unity throughout this paper. 

This matrix cannot reproduce the realistic structure of CKM matrix and quark masses although 
it seems roughly consistent with them. Since this matrix is singular, we cannot have a definite 
unitary matrix diagonalizing Y -Y^ . In addition, the resultant CKM mixing matrix does not have 
any mixing between the third generation and the other ones. Besides, only one of three eigenvalues 
is non-zero for each sector, which means that the model cannot explain the masses of quarks and 
leptons except for the third generation. As for this problem, several attempts have been done for 
generating definite non-zero mixings. For example, the authors in Ref. [20, 21] argue on the basis 
of supersymmetric composite models constructed in the framework of intersecting brane model. 
The authors in Ref. [26] study the CKM mixing matrix in the model with 3 supersymmetric Higgs 
doublets, where all the matter fields arise from one of the tori. 

Let us assume here that the metric of the first torus depends on a coordinate on the third torus 
xg. Then, the length of strings projected on the first torus become dependent on the coordinate 
xg, and so does the area swept by the strings. In this case, Yukawa couplings cannot be expressed 
in the simple form (3.15). In a rough approximation, we take the effect of the xg-dependence into 
consideration by modifying ^4^^ 

Af^Af{l-ei), (3.17) 

where is a small quantity generated by the xg-dependence. The is, in principle, calcula- 
ble in string theory, but we treat them simply as input parameters. In this estimation, Yukawa 
couplings are no longer of the form (3.15), and thus Yukawa matrices would generically become 
non-degenerate. Surprisingly, this slight modification gives the realistic hierarchical structure of 
fermion masses and a CKM mixing matrix. In this paper, we attribute the xg-dependence of the 
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metric to a warp effect of the background geometry.^ Before calculating quark masses and CKM 
mixings, we here make a comment on the possibility of the emergence of warping in the intersecting 
brane model introduced in the previous section. 

The warped geometry was first proposed in Ref. [30], where the authors began with a five 
dimensional model M'^ x S^/Z2, with cosmological constants in the 5D bulk space and on the 
4D boundaries, and found a solution of the Einstein equation which shows that the 4D metric 
depends on the extra dimensional coordinate 0: ds"^ = e~'^^^''^'^^rinydx^dx'^ + r'^dcf)^ , where A: is a 
five dimensional scale and is the compactification radius. 

It must be noted that the warping comes from localization of the energy density. In our setup, 
the internal manifold is T^x 

ji2 x2^2 g^j^^ orientifold planes are located at a;^ = 0, vr; = 0, vr; x^ = 
0, vr, and D-brane stacks are wrapped on three cycles on the tori. For simplifying the situation, let 
us smear the D-branes and orientifold planes on the first and second torus, not smearing them on 
the third torus, and assume that the D-brane stacks projected on the third torus are all parallel 
to the orientifolds. Looking at this configuration neglecting all the extra dimensions except for x^^ 
one could regard it as the 5 dimensional model compactified on /Z2 with two negative tension 
branes and several positive tension branes. Because of the localization of the positive energy branes 
and the negative energy branes, it seems valid to assume the emergence of a warp factor depending 
on x% though it might be small. 

In the present paper, we do not attempt to solve the metric background, but it is interesting to 
note that the authors in Ref. [32] investigated a similar setup. They consider 10 dimensional Type 
IIA string theory compactified on IZ2 with parallel 32 D8-branes and 08-branes. Assuming that 
2n D8-branes are located on the top of the orientifold plane at x% = 0, and 32 — 2n D8-branes on 
the orientifold plane at xg = ttR, where R is the radius of S^/Z2, they found the solution 

ds^ = H-l^^l-dt^ + {dx^f] + Hl^^^{dx^f, (3.18) 

where i^DS = 1 ~ ^Dsl^^^j; ^D8 = ^'^2^1^" ' 3^ string coupling constant, and Ig = Voi' is the 

string length. If we compactify this configuration on tori and T-dualize it twice, it becomes close 
to our model. Although the procedure of torus compactification might have an influence on the 
warp factor in (3.18), it is not so audacious to expect in our setup some kind of warping along one 
of the internal coordinates. 

Now let us move to the evaluation of Yukawa matrices and calculate the CKM mixing matrix. 

^The dependence might be explained by other mechanism. For example, in the context of M-theory compactified 
on /Z2 X Ke, where Kq is a Calabi-Yau manifold, the authors in Ref. [31] solved the equation of motion of IID 
supergravity by using a strong coupling expansion, and found that the volume of the Calabi-Yau manifold can linearly 
depends on the coordinate of jZ^- 

*The authors in Ref. [33] suggest that when one would like to realize a five dimensional Randall-Sundrum type 
solution from higher dimensional theory such as supergravity, M-theory or F-theory, it would be more practical to 
first solve the Einstein equations and obtain a Randall-Sundrum type solution in higher dimensional spacetime and 
then perform some dimensional reduction than to do reversely. 
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Turning on the warp effect, we have Yukawa matrices for the up-sector: 



A% = Af + A^l-xci). 



(3.19) 
(3.20) 



In this expression, we put in the form = xq, where a; ^ 1 and q is of order 1 since the order 
of ei is thought to be of the same order. In order to calculate the masses and the CKM matrix on 
the basis of series expansion, besides x, we define two more small quantities y and z as 



-kA2 



-kAi 



Z = 



{Ai > A2> A3). 



(3.21) 



^ Q—kAs ' Q—kA2 ' 

Solving the eigensystem of Yukawa matrices by series expansion with respect to x, y and z, we 
finally obtain the leading part of up- and down- quark masses 
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and the CKM mixing matrix 
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(3.23) 
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(3.30) 
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i=l 



(3.34) 



If we ignore coefficients and assume x, y and z of order e, we can reproduce the hierarchical structure 
of the realistic CKM matrix: 

9 
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(3.35) 
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In this naive estimation, the masses of the first and the second generations, mi ^ and m2 ^ xy, 
turn out of the same order e^, and one might regard it as an undesirable result. However, noting that 
X comes from a small warp factor, and that it should be smaller than y and z, the mass difference 
is expected to come from a small value of this ratio x/y and, of course, from the coefficients we 
have so far neglected. 

4 Discussion 

In this paper, we propose a method for reproducing the realistic hierarchical structure of the 
observed quark masses and CKM mixing matrix on the basis of a specific intersecting brane model 
proposed in Ref. [9]. We newly introduced a small warp factor on tori, which systematically 
changes the elements of hierarchically suppressed but singular Yukawa matrices. The modification 
gives two hierarchically suppressed quark masses and that of order 1 in addition to a definite CKM 
mixing matrix. Taking the order of several parameters appropriately, we reproduce the realistic 
hierarchical structure of quark masses and a CKM mixing matrix. It is worth noting that both 
large mass hierarchy and small mixing angles are originated from a warp factor on the internal 
manifold. 
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